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Abstr ict. We present basic laboratory infrared data on et al. 1995|) and recent g round-based observations with 



a large number of SiC particulate samples, which should 
be of great value for the interpretation of the 11.3 /im 
feature observed in the spectra of carbon-rich stars. The 
laboratory spectra show a wide variety of the SiC phonon 
features in the 10-13 /im wavelength range, both in peak 
wavelength and band shape. The main parameters de- 
termining the band profile are morphological factors as 
grain size and shape and, in many cases, impurities in 
the material. We discovered the interesting fact that free 
charge carriers, generated e.g. by nitrogen doping, are a 
very common characteristics of many SiC particle sam- 
ples. These free charge carriers produce very strong plas- 
mon absorption in the near and middle infrared, which 
may also heavily influence the 10-13 /jm feature profile 
via plasmon-phonon coupling. 

Wei also found that there is no systematic dependence 



of the I band profile on the crystal type (a vs. (i SiC). This 
is proven both experimentally and by theoretical calcula- 
tions based on a study of the SiC phonon frequencies. 
Further, we give optical constants of amorphous SiC. We 



UKIRT flSpcck ct al. 19971) . The UKIRT data showed that 
the feature cannot only appear in emission, but that there 
are a few sources with an absorption feature. 

An ultimate proof for the formation of SiC grains in 
carbon-rich stellar atmospheres was the detection of iso- 
topically anomalous SiC grains in primitive meteorites 
(Bernatowicz et al. 1987). The s-process isotopic signa- 
ture identified in the majority of the meteoritic SiC grains 
can be produced in the atmospheres of thermally pulsing 
low-mass AGB stars ( Gallino et al. 1994 ). Most of these 
presolar SiC grains seem to be large with equivalent spher- 
ical diameters between 0.3 and 3/im ( Amari et al. 1994 ), 
but also very large grains up to 26 /Ltm are present. Such 
large grains can probably only be formed during extreme 
mass loss phases or in clumps of high density ( Bernatow- 
1996). So far there is no spectroscopic evidence 



et al. 



for SiC grains in the diffuse interstellar medium. This fact 
is not yet satisfactorily explained although core-mantle, 
inhomogeneous, or very large grains may offer a solution 



discuss, the implications of the new laboratory results for (fosenkopf et al. 1992;; [Kozasa et al. 1996|; |Andersen et al 



the int erpretation of the spectra of carbon stars . 



1998) 
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1. Introduction 

Thermodynamical equilibri um calculation s performed by 
Friedemann (1969a, b) and Gilman (1969) suggested that 
silicon carbide (SiC) particles can form in the mass out- 
flows from carbon stars. First empirical evidence for the 
presence of such particles came from infrared spectroscopy 



The 11.3 /im feature in carbon star spectra displays 
some interesting variability which may offer the possi- 
bility to characterize more closely the composition and 



structure of cosmic SiC grains. Littlc-Marcnin (1986) an- 
alyzed 176 spectra of bright infrared sources in the IRAS 
LRS database and found for 4% of the carbon stars a fea- 
ture shifted to w 11.6 /im. The majority of the features 
peaked at 11.15 ± 0.10 /im with a FWHM of 1.6 ± 0.15 
/im. The statistical analysis of spectra of class 4n revealed 
some more variety with only part of the spectra peaking 
at 11.3 /im and a second peak at 11.7/im appearing when 



perform ed by Hackwell (1972) and Trrjffeis fe Cuhen 



the ratio of feature to continuum decreases (Baron et al 



(1974) The latter anLhuis iiiLeipieLed the deLecLed iul- 



1987). A very thorough study of 718 LRS spectra of carbon 



atively broad 11.3 /im feature as the emission of small SiC 
particles between the transverse and longitudinal optical 
phonon frequencies. More empirical material was accumu- 
lated by the IRAS mission ( Little- Marenin T986] ICocbcj 
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stars based on autoclassification methods was performed 
by Goebel et al. (1995)| who also found a weakening, broad- 
ening and shift (to 11.9 /im) for classes of objects with a 
low conti nuum temperature. The ground-based data ob- 
tained by Speck et al. (1997) also show in about one third 
of the features a peak at wavelengths larger than 11.4/im, 
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but without a strong correlation to the blackbody tem- 
peratures. 

A definite explanation of these spectral differences in 
terms of an identification of a well-defined silicon carbide 
material as carrier of the 11+ fim feature would be very 
valuable, because it could provide important information 
concerning grain condensation mechanisms and the condi- 
tions within circumstellar grain- forming regions. Unfortu- 
nately, this task is complicated by (i) the presence of only 
one relatively broad band, (ii) some uncertainty in the 
location of the continuum due to nearby molecular bands 
( Hron ct al. 1998| ) (iii) possible contribution from carbona- 
ceous material to the emission, (iv) the extreme sensitivity 
of the SiC band profile to particle size and shape effects, 
(v) the appearance of many different SiC crystal types 
(polytypes), (vi) the possibility of disordered structures, 
and (vii) the influence of material impurities on the band 
profile. 

Three different explanations of the observed band vari- 
ety have been given so far. Baron ct al. (1987) and |Goebc] 
et al. ( )-995)| found a correlation of the apparent band shift 
to the appearance of a feature at about 8.5 /xm as well as to 
an increase in the FIR flux and attributed these combined 
spectral changes to an increasing contribution by hydro- 
genated amorphous carbon (a:C-H). |Speck et al. (199*7) 



explained the different band profiles by different crystal 
types of the SiC grains and tried to fit the UKIRT spectra 
to laboratory spectra (see below). |Papoular ct al. (1998) 



showed that also different grain shapes can produce the 
observed variation in band shape. 

Considering this complexity of the identification prob- 
lem (see also Sandford 1996), detailed laboratory investi- 
gations on the infrared spectrum of SiC particles unfor- 
tunately are scarce both in the physical/chemical and the 
astronomical literature. Most popular among astronomers 
are the studies by Stephens (1980) on laser-produced f3- 
SiC condensates, 



Friedemann et al. (1981) 



cially available a-SiC, and Borghesi et al. (1985) 



on commcr- 
on com- 



mercially produced a- and /3-SiC. More recently, Papoular 



et al. (1998) investigated two samples of /3-SiC powders, 
one produced by laser pyrolysis and one which was com- 
mercially available. In the earlier papers, influential fac- 
tors beyond the polytype remained mostly unexplored. Es- 
pecially, the decisive influence of grain shape and size as 
well as of the surrounding medium (ma trix) on the band 



profile, which was already discussed by Treffers fc Cohen 
(1974)[ and later investigated in detail by Bohrcn fe Huff- 



man (1 983) , has not been treated correctly or even ignored 



by the majority of the authors. Papoular ct al. (1998) 



gently brought this topic back to mind and presented an 
exact method for the correction of matrix effects. 

As already pointed out, the available experimental 
studies with a total of about 15 samples are by no means 
systematic. Moreover, the lack of good laboratory results 
led to some confusion concerning the influence of the crys- 
tal type on the band profile. The goal of our paper is to 



present experimental data on a-SiC, /3-SiC, and amor- 
phous SiC and to investigate experimentally size, shape, 
and matrix effects (Sect. 3 and 4). To give the discussion 
on the optical properties of different SiC polytypes a more 
exact foundation, we summarize in Sect. 2 basic facts on 
SiC bulk crystals and thin films and calculate theoreti- 
cal spectra of small SiC grains for different polytypes. In 
Sect. 5 we discuss in which way the new data can serve 
the observers for an evaluation of their data in terms of 
the SiC hypothesis. 

2. Basic facts on structural and infrared optical 
properties of SiC 

2.1. Crystal structure of SiC 

The crystal structure of SiC shows pronounced polytypism 
which means that there exist a number of possible crystal 
types differing in only one spatial direction (see Fig. [l]) . 
The basic units from which all polytypes are built are 
Si-C bilayers with a three-fold symmetry axis, in which 
the Si- and C-atoms are closely packed. These bilayers are 
stacked on top of each other while they are laterally shifted 
by 1/V3 of the Si-Si or C-C atomic distance in the layer 
either in the [1100]- or in the opposite direction. Hence, 
each Si atom is tetrahedrally surrounded by four C atoms 
and vice versa. If all shifts occur in the same direction, 
then an identical position of the bilayer in the projection 
along the hexagonal axis is reached after three stacking 
steps. The resulting structure is of cubic symmetry and 
because of the three-step stacking period this polytype is 



called 3C (C for cubic, Ramsdell 1947). Another name for 



this polytype, which is the only cubic one, is the often- 
used term /3-SiC. The environment of each Si-C bilayer 
in the stack, which is determined by shifts in the same 
direction at each step, is also called cubic. 

The other extreme is obtained, when the bilayers are 
shifted alternatingly in opposite directions such that, in 
projection with the hexagonal axis, every other layer has 
the same position. The lattice is then of hexagonal type, 
and because of the two-step period the polytype is called 
2H. The environment of each Si-C bilayer which has been 
produced by the alternating shiftes is also called hexago- 
nal. 

All other polytypes are built up by a characteristic 
sequence of cubic and hexagonal Si-C bilayers, for which 
the 3C and 2H polytypes represent the limiting cases. All 
polytypes except 3C are uniaxial crystals (optical axis = 
c-axis) and belong either to the hexagonal or to the rhom- 
bohedric system. The most abundant polytypes besides 
3C and 2H are the hexagonal types 6H and 4H and the 
rhombohedric 15R. For historical reasons any non-cubic 
polytype or mixture of polytypes is also designated by the 
term a-SiC. 

The ratio of the numbers of hexagonal to cubic bilayers 
is called hexagonality and is a very useful scaling param- 
eter. Several properties of the polytypes change with this 
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Fig. 1. Primitive hexagonal unit cells of the most simple 
SiC polytypes. Si atoms are represented by open circles, 
C atoms by filled circles. Bilayers of the three possible 
positions in projection with the c-axis are labeled by the 
letters A, B, and C. The Si-C bonds in the (1120) plane 
indicating the relative shifts of the bilayers are represented 
by heavy solid lines. The figure has been reproduced from 
Kackell et al. (1994) . 



parameter, e.g. str uctural propert i es such as the lattic e 
constants (see, e.g., Madelung 1982 ; Bcchstcdt et al. 1997 ), 



the energy o f the indirect band gap ( Feitknecht 1971 



Choyke 1990), and certain phonon frequencies (Nakashima 
et al. 1 987| ; Hofmann et al. 1994 ) , which will be discussed 



below. 



2.2. Lattice vibrations of bulk SiC and thin films 

As in any other ionic crystal, the interaction of SiC with 
infrared light is clearly dominated by the fundamental vi- 
brations of the Si- and C-sublattices against each other. 
The excitation of these phonons produces strong fea- 
tures in the wavelength range 10-13 /im, observable e.g. 
in reflectance measurements as a nearly 100% reflectiv- 
ity (so-called reststrahl band). The absorption coefficient 
(depending on the crystal quality) can be up to several 
10 5 cm -1 . Therefore, transmission measurements in this 
spectral region are only possible by using samples with 
effecth p thicknesses below some 100 nm. 

In general, the trcqucncics ot the involved longitudinal 
(LO-) and transverse (TO-) optical phonon modes depend 
on their polarization direction in relation to the c-axis. 
One reason for this behaviour - in simplified words - is 
the slight stretching of the basic tetrahedra in c-dircction 
causing a smaller force constant in hexagonal double lay- 



limiting case are the modes polarized parallel to the c-axis 
(Ai-modes), which have a smaller resonance frequency 
than the Ei modes. For the transverse modes, the dif- 
ference is nearly proportional to the fraction of hexago- 
nal double layers (the hexagonality parameter h) in the 
polytype. This difference was experimentally found to be 
gj(E 1T ) - m(A 1T ) = 29 .4 cm" 1 - h ( |Feldman et al. 1968b ; 
jNakashima et al. 1987 ) and theoretically calculated with 
lu(Ei T ) - lu{A 1t ) = 29.2 cm" 1 - h (|Hofmann et al. 1994|). 



Averaged values of published data of the mode fre- 
quencies of some polytypes, determined by Raman and 
infrared measurements, are given in Table 0. For compari- 
son, theoretical values derived from the bond charge model 
(Hofmann et al. 1994) are also given. For the 3C polytype, 
the maximum deviations of the single published frequen- 
cies from the average values in Table [j] are 1.9 (TO) and 
1.2 cm -1 (LO). The biggest deviations for the Eit, Ait, 
Ei/,, an d Ail modes of the 6H polytype are 1.0, 0.9, 1.3, 
and 1.7 cm , respectively. 

~~ for the TO- and LO 



The data by Spitzer et al. (1959a 



mode frequencies of the 3C polytype, which have been of- 
ten used in astronomical papers, seem to deviate system- 
atically from the mean values (-1.9 cm -1 and -2.0 cm -1 , 
respectively). For the Eit and Ai£ mode frequencies of 
the 6H polytype ( |5pitzer et al. 1959b| ), the deviations are 
even bigger (-3.1 and -4.6 cm -1 , respectively). These de- 
viations probably indicate limited accuracy either of the 
measurements or of the fits in this work compared to newer 
ones. Therefore, we exclude these results from the the av- 
eraging. Including them would result in shifts of the mean 
values by an amount of up to 0.4 cm- 1 . 

Depending on the complexity of the polytype struc- 
ture, there are a number of weak modes in addition to the 
strong modes discussed so far. With increasing unit cell 
length and, therefore, increasing number of atoms in the 
unit cell, the number of branches in the Brillouin zone in- 
creases. Because these branches can be folded back to the 
phonon dispersion curve of the basic 3C-SiC polytype with 
the largest Brillouin cell ( Patrick 1968| ), these polytype- 



specific highly anisotropic modes are called zone-folded 
modes. They can be observed by Raman spectroscopy (e.g. 



Feldman et al. 1968a; 


Nakashima & Harima 1997), in- 


frared reflectance (e.g. 


Tin & Rashevskaya 1973; Pcnsl & 



Helbig 1990) and transmittance measurements with very 



ers compared to the cubic ones (Nakashima et al. 1987). 
One limiting case are the modes polarized parallel to the 
basal plane (perpendicular to the c-axis, Ei modes), which 



are veriy similar fur Lhc different pul y types. Their frequen- 
ci es differ nut more Lhan 3-4 cm- 1 (Feldman cL al. 19081) 



thick samples (e.g. Dubrovskii & Radovanova 1973). 

The peak absorption coefficient caused by these folded 
modes is approximately 10 2 cm -1 . Because this is three 
orders of magnitude lower than the value for the strong 
bands, the weak modes are practically not observable in 
astronomical spectra. The same is true for the combina- 
tion modes of lattice vibrations which cause additional 
weak absorption effects e.g. in the spectral re gion 1250- 
1700 cm -1 (two-phonon combination modes, Patrick & 



Nakashima et al. 1987j |Hofmann et al. 1994| ). The other 



Choyke 1961) 
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Table 1. Hexagonality parameter h, resonance frequencies of the strong phonon modes, and high frequency dielectric 
constants of some SiC polytypes. Eit and En, denote the frequencies of the TO- and LO-modes in the basal plane of 
the SiC lattice, respectively, whereas Ait and Ail denote those parallel to the principal axes. 
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2.3. Modelling the infrared optical properties of SiC 

SiC can be regarded as a textbook example concerning the 
modelling of reflectance spectra of bulk and thin film sam- 
ples (e.g. Bohrcn & Huffman 1983; Mchnichuk & Pascch 



nik 19p[ |Engclbrccht fc Hclbig 1993|) and of trans mit- 
tance spectra of thin films (e.g. Spitzer et al. 1959a) by 
a Lorentz-oscillator model of the dielectric function. In 
many cases it is sufficient to regard only the strong modes 
so that the dielectric function parallel and perpendicular 
to the optical axis can be parametrized by a simple one- 
phonon Lorentz term: 



with 



-JCO 



u jTO 



J jTO 



) (i=IU)- 



(i) 



(2) 



Here the quantities uj±to ( w _llo) andw||TO ( w ||lo) are 



^it anu ait (A 1L 

respectively, tire plasma, frequencies u jP represent Lhe us 

cillator strengths and the jj describe the phonon damp- 



the visible and NIR spectral range. Some average exper- 
imental and theoretical values of e JOO from the literature 
are given in Table |l|. The scattering of the experimental 
data in the literature is relatively large (up to ±0.2) but 
the agr eement of the avera ge values with the theoretical 
data by Chen et al. (1994) is satisfactory. 

The damping constant 7 (inverse proportional to the 
phonon life time) is an "ad hoc" introduced parameter, 
which in a perfect crystal reflects the anharmonicity of 
the potential curve. On the one hand, the phonon lifetime 
is limited by coupling processes, opening energy dissipa- 
tion channels. On the other hand, phonons are scattered at 
crystalline defects. Therefore, the damping constant can 
be taken as a measure of crystallinity. High-quality SiC is 
characterized by damping constants of 1-3 cm -1 . These 
values are reached by bulk crystals (Peiter, unpublished) 
as well as by micron-thick CV D-grown layers pro duced at 
relatively high temperatures ( Zorba et al. 1995 ). Single- 
crystalline layers grown by solid-source MBE at lowe r tern 
peratures show damping constants down to 5 cm -1 (Pfcn 



the frc ^iicios of Llic E 1T (Eg and A 1T (A 1£ ) modes, nighaus ct g ^ Polycrystalline layers pro duceTfo7 



ing (liim width) of the corresponding mode. The quantity 



instance by sputtering or laser ablation require damping 
constants of several 10 cm -1 in the modelling (Hobert 



e JOO is 1 Lhe high-frequency dielectric constant caused by 

the system of the valence electrons, which is measured in 



et al. 199E) 
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Fig. 2. Calculated absorption cross section per unit volume C/V for small SiC grains of different shapes and polytypes. 
Profiles for spherical grains are given by solid lines, the CDE according to Bohren & Huffman (f983) by dotted lines, 
and a CDE with higher probability for less elongated grains (Ossenkopf et al. 1992) by dashed lines. 



2.4- Influences of real structure and impurities 

In the preceding sections more or less perfect SiC crystals 
have been considered. However, for SiC grains grown by 
kinetically controlled reactions in the outflow of carbon 
stars it is important to take the influences of structural 
disorder, non-stoichiometry, and impurities on the infrared 
absorption into account. 

Great experience on these effects is already available 
from thin film studies. The influence, e.g., of disorder up to 
complete amorphousness, has been investigated with films 
produced either by implantation of carbon ions into silicon 



wafers (e.g. Kimura et al. 1981; Serre et al. 1996) or by 



amorphization of crystalline SiC by ion irradiation ( Zorba 
et al. 1 996| ; Musumcci et al. 1997). The phonon band pro- 
files of these films in some cases still have been described 
by the Lorentz oscillator model with a very large damping 
consta nt of about 200 cm' 1 ( Zorba et al. 1996|). However, 



Musumcci et al. 1997 ). The modelling of these spectra re- 
quires a Gaussian distribution of oscillators. Additionally, 
the centre of the oscillator frequencies shifts to smaller 
frequencies by several 10 cm -1 . The reasons for these ef- 
fects are the variation and the (on average) increase of 
the bond lengths in the amorphous solid compared to the 
crystalline one. Additionally, the break-down of the k«0 
selection rule causes contributions from the whole optical 
branches to the absorption spectrum. Both the broaden- 
ing and the shift of the absorption feature can be avoided 
or inverted by high temperature during or after the sam- 
ple preparation, which supports recrystalization (see, e.g., 



Durupt et al. 1982) 



Non-stoichiometry, especially carbon deficiency, leads 
to the occurrence of separated SiC grains in some kind 
of matrix. The phonon band observed at such samples is 



shifted to h igher frequencies (Yoshii et al. 1991; Lundquist 



very o ten a change of the band profile from Lorcntzian 
to Gaussian shape is observed (e.g. Basa fc Smith 199C ; 



et al. 1995 ), which has been attributed to surfac e excita- 



tions of the SiC grains ( Bean fc Newman 1971 ). In con 
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trast to this, carbon excess leads to additional continuum 
absorption caused by the semi-metallic character of the 
graphitic excess carbon. 

Impurities and defects incorporated into the SiC crys- 
tal prod uce besides a slight increase of the p honon damp- 
ing (e.g. [Tin et al. 1972 ; Harima et al. 1995) the following 
effects: a) localized vibra tional modes ( Newman 1990 ; En 
gelhrec ht fe Helbig 1994), b) photoionization bands of im- 



purity atoms (e.g. Dubrovskii fc Radovanova 1971 ; Purt 



seladzc fc Khavtasi 1972| ), and c) free charge carrier ab- 
sorption. The resulting absorption coefficient from these 
impurity-induced effects is in bulk material and films usu- 
ally two or more orders of magnitude smaller than the 
contributions from the strong phonons. The free charge 
carrier absorption, dominant at longer wavelength and 
higher temperatures, depends on the degree of doping 
( Kulakowskii et al. 1975| ) and may become so strong that 
it can be observed even in the mid-IR reflectance (carrier 
density >4-10 17 cm -3 ) ( |lmai 1966 ; Mcl'mchuk & Pascch 
nik 19^). The longitudinal optical phonon and the free 
charge carrier collective oscillation (plasmon) generally 
form coupled modes which occur in two branches (L + 
and L_) at the high and low frequency sides of the rest- 
strahl region flHolm ct al. 1986| |Sasaki et al. 1989| ). The 
modelling of such spectra requires an additional oscillator 
term representing the response of the free charge carrier 
gas to electromagnetic excitation (so-called Drude term) 
in Eq. (|l|) . This Drude term has a similar structure as the 
Lorentz term, but with the resonance frequency being zero 
and the plasma frequency being proportional to the square 
root of the charge carrier density. 



2.5. T heoretical absorption spectra of very small spherical 
and ellipsoidal tiiU grams of different polytypes 

In contrast to bulk material, the vibrational absorption 
bands of small grains are determined by surface modes 
which occur between the LO- and TO-frequencies of 
strong modes, the exact position dependi ng on size, shape, 
and the medium surrounding the grains dBohrcn fc Huff] 



man 1983). Within the limit of particles very small com- 
pared to the wavelength (quasistatic or Rayleigh limit), 
there exist simple and exact formulae for the calculation 
of absorption spectra of ellipsoidal grains from the dielec- 
tric function. An ensemble of such (identically shaped) 
grains embedded in random orientation in a matrix with 
the dielectric constant e rn will cause three sharp spec- 
tral features (surface modes) at the resonance frequencies 
w i = ^to + Liujl/(e m + Li(eoo - e m )), (i=l,2,3) where 
the geometrical factors Lj (Lj >0, ^Lf=l) characterize 

i 

the lengths of the principal axes of the ellipsoids (L;=l/3 
for spheres, 0, 1/2, 1/2 for needles, 0, 0, 1 for discs, see also 



Papoular et al. 1998). Unfortunately, for arbitrary, espe- 



cially sharp-edged shapes, theory is still far from the abil- 
ity to determine the absorption profile caused by surface 
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Fig. 3. Calculated absorption cross section per unit vol- 
ume for small a-SiC grains (6H polytype) taking into ac- 
count the weak zone-folded modes. The lines designate 
different grain shape distributions (compare Fig. ||) . For 
clarity only the polarization direction parallel to the c- 
axis is shown. 



modes of grains. A widely used possibility for providing 
a rough estimate for the absorption spectrum of a col- 
lection of real particles is the continuous distributions of 



ellipsoids (CDE) introduced by Bohren fc Huffman (1983 



Besides this (rather extreme) uniform distribution of all 
possible ellipsoidal shapes we also use a variation of the 
CDE which considers a weighted distribution with max- 
imum probability for spheres (m-CDE, pssenkopf et al 



1992) 



From the dielectric function Eq. ([!]) with the resonance 
frequencies and high-frequency dielectric constants given 
in Table [l] we have calculated theoretical infrared absorp- 
tion spectra for small SiC spheres as well as for the two 
CDE approximations, all in a vacuum environment. The 
results of these calculations for the polytypes 3C (/3-SiC), 
2H, 4H, and 6H are given in Fig.||. A damping constant 
7=2 cm -1 which is characteristic of a good-quality crys- 
tal (see Sect. 2.3) was used. The orientations of the crystal 
axes with respect to the polarizations of the incoming ra- 
diation have been averaged by considering one third of the 
axes being oriented parallel and two thirds normal to the 
polarization direction. 

The calculated band profiles show the expected strik- 
ing dependence on the grain shape distribution. For 
spherical grains very sharp resonances occur, whereas 
for the CDEs the wide distribution of shapes re- 
sults in broad bands between the LO and TO fre- 
quencies. For the anisotropic modifications (the a-SiC 
polytypes 2H, 4H, and 6H), the two principle orien- 
tations produce their own resonances in the sphere- 
spectra. The positions of the resonances are 10.71 /zm 
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Fig. 4. Calculated absorption cross section per unit vol- 
ume for small a- and /3-SiC grains using a higher damping 
constant of 7=10 cm -1 than in Fig.|] (2 cm -1 ). 



for the 3C-, 10.73 /im and 10.77 fim for the 6H-, 
10.75 /im and 10.80 /im for the 4H-, and 10.68 /im and 
10.78 /im for the 2H-modifications, respectively. For com- 
parison, the resonance positions using the data by Spitzer 
et al. (1959a,b) would be 10.72 /tm for 3C-SiC, and 
10.73 /im and 10.84 /im for 6H-SiC. 

The CDE spectra for the polytypes essentially differ 
in a shoulder at the long-wavelength side of the profile, 
which is also caused by the anisotropy. Since the TO- 
frequencies in the two principal orientations deviate from 
each other, the band profile extends to different spectral 
positions which results in a shoulder after the orientations 
have been averaged. The difference of the E%t and Ait 
frequencies and consequently the broadness of the shoul- 
der is correlated to the hexagonality in the sense that for 
the 6H modification the shoulder is only 0.14 /im broad 
but for the 100% hexagonal 2H modification its width is 
0.48 /im. In the m-CDE the shoulder is not present since 
the extreme shapes producing it have about zero proba- 
bility. 

So far we have only considered the strong phonon 
modes. To check if the weak "zone-folded" modes could 
deliver a criterion useful for distinguishing the polytypes 
in small-particle spectra, we show in Fig.|^ as an exam- 
ple a calculated absorption spectrum of the 6H poly- 
type, which includes the weak modes. These have been 
represented by two additional oscillators with the TO- 
and plasma-frequencies W2to=883.8 cm -1 , W2P=50 cm -1 , 
cj3TO=888.6 cm -1 , uj^p=A5 cm -1 and with 7=2 cm" 1 
for both oscillators. The effect of the weak modes is 
the appearance of two very weak features at their TO- 
frequencies, interestingly as dips in the CDE spectra. Note 
that after averaging with the two times stronger contribu- 
tion of the perpendicular orientation these features will be 



even weaker and hardly detectable in astronomical spec- 
tra. No influences on the resonance positions and the lim- 
iting frequencies of the band have been found. 

Fig. |4] shows calculations using a higher damping con- 
stant of 7=10 cm -1 which would characterize less struc- 
turally perfect crystals but material still far from amor- 
phousness. The calculations have been done for the most 
abundant a-SiC polytype (6H) and for /3-SiC. They show 
that for such less perfect crystals the anisotropy effects in 
the a-SiC spectrum are nearly smeared out. The double 
feature of the spheres merged into one band with a peak 
position very close to that of the /3-SiC spheres. The shoul- 
der in the CDE spectrum at u>to changed into a smooth 
slope which is only a little less steep than for the /3-SiC. 

Therefore, the result of the theoretical considerations 
presented in this section is that the small differences be- 
tween the phonon frequencies of the SiC polytypes lead 
to detectable features in small-grain spectra only if the 
crystal structure of the grains is very perfect. Even in this 
case it is required that either the grains have a very defi- 
nite shape (e.g. spherical) in order to produce sharp res- 
onances, or that there is a sufficient amount of very elon- 
gated structures to produce enough absorption at the TO 
frequency so that the shoulder caused by the anisotropy 
of this frequency can be seen. The latter seems not to- 
tally unrealistic. However, we would like to stress that in 
any case the differences will be more marginally and that 
(in agreement with Papoular et al. 1998| ) the general ap- 
pearance of the feature as well as the peak position only 
depend on the grain shape. 



3. Transmission measurements on SiC particles 

3.1. Materials and measurements 

We studied 16 different SiC powders which are partly of 
commercial origin and partly laboratory products (see Ta- 
ble [2]) . The names we assigned to these samples serve only 
for an easier distinction. For some of the commercial pow- 
ders, information about the production route and the pu- 
rity is very poor. Catalogue values of purity take in some 
cases only metal impurities into account and not the ma- 
jor pollutants: free C, Si and Si02- 

Generally, the samples fall into two groups. The first 
one is the commonly used SiC for grinding application, 
which is produced from quartz and coal in an electric fur- 
nace at very high temperatures (about 2500°C, Acheson 
process) . These powders are often large in grain size (Pies- 
teritz, Duisburg) but ESK, Lonza and probably AJ-a also 
belong to this group. All of these high-temperature prod- 
ucts consist of a-SiC (dominated by the 6H polytype). 

The second group is produced from gaseous or liquid 
precursors in some kind of pyrolysis which may be induced 
by a plasma discharge, laser irradiation or simply by heat- 
ing. This kind of process is capable of producing very fine 
(submicron) powders and allows easy mixing with other 
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Table 2. Properties of the investigated SiC samples. Polytype according to X-ray diffraction, size and shape from 
optical or transmission electron microscopy. Regular shape means observation of hexagonal/triangular (Dresden, Plas- 
machem) and spherical (Jena-Ar) grain contours. Purity is given according to elemental analysis (partially provided 
by producer); question marks indicate that only metal impurities have been considered. Major impurities are given 
in parentheses. Production method Acheson indicates high-temperature products, GPP indicates gas-phase pyrolysis 
(see Sect. 3.1.) 



product name 


polytype 


max. size 

(H 


shape 


purity 
(wt.%) 


prod.meth. 


AJ-a 


a 


(6H) 


2 


irr. 


99.8? 


probably Acheson 


ESK 


a 


(6H) 


2 




98 


Acheson 


Lonza UF-15 


a 


(6H) 


1 






probably Acheson 


Duisburg 


a 


(6H) 


40 


irr. 




probably Acheson 


Piesteritz (green) 


a 


(6H) 


150 


irr. 


98.8 


Acheson 


Piesteritz (black) 


a 


(6H,4H) 


150 


irr. 


96.4 (C) 


Acheson 


AJ-/3 


P 




1 


irr. 


96.5 




AJ whiskers 


P 




2x60 


wh. 


99.5? 




Plasmachem 


a 


(6H) 


0.5 


reg. 


95.2 (Si0 2 ,Si 3 N 4 ) 


plasma 


Starck B20 


P 




0.6 


irr. 


95 




Japan 


P 




0.5 


irr. 




GPP 


Dresden 


a 


(6H) 


0.3 


reg. 


83.6 (Si 3 N 4 ,Si) 


GPP, polysilazane 


Riga 


P 








86 (Si 3 N 4 ,Si) 


plasma 


Berlin 


P 




0.5 


irr. 


85 (C) 


plasma, SiCU 


Jena 


P 






reg. 




laser pyr., SiH 4 ,C2H2 


Jena-N 


P 






reg. 




laser pyr., SiH 4 ,C 2 H 2 ,NH 3 



components such as S13N4. Therefore, a huge variety of 
such processes and products has already been studied in 
the literature. The polytype of these grains is usually 3C 
(/3-SiC), but depending on the process temperature also 
a-modifications (higher T, see Dresden, Plasmachem) or 
even amorphous components may occur. The purity is 
mostly not worse than that of the high-temperature prod- 
ucts (if not otherwise desired as the S13N4 component in 
Dresden and Riga) and the major impurities are the same 
(free C, Si, oxide). For the sake of comparison we also have 
one product with a big carbon excess (Berlin). 

For astrophysically motivated studies, SiC samples 
from pyrolysis processes are certainly preferabl e to the 
high-temperature products. 



Cherchncff (1997) 



investi- 
gated SiC particle formation in stellar atmospheres and 
showed that hydrogen must play an important role in the 
formation chemistry. In hydrogen-deficient environments 
as R CrB and Wolf-Rayet star atmospheres, SiC emission 
is not observed. The formerly proposed formation mecha- 
nism of SiC in carbon stars by reaction of SiC2 and Si was 
shown by Willacy fc Cherchncff (1998)| not to be efficient 
enough to produce the observed quantities of SiC dust. 
If this is true, the pyrolysis processes which in nearly all 
cases start from hydrogen compounds of C and Si, will cor- 
respond much better to the astrophysical grain formation 
process. 

For all powders we determined the polytype(s) by X- 
ray diffraction as well as grain shape and size by opti- 



cal or transmission electron microscopy (Table §. For the 
Piesteritz products which have already been studied by 
Fricdcmann ct al. (1981), the X-ray analysis revealed a 



simpler composition than that supplied by the producer 



which was published by Fricdcmann ct al. (1981 



Infrared spectra in the wavelength range 2-25 fim have 
been recorded with a Bruker 113v FTIR spectrometer on 
the SiC powders embedded in KBr pellets with a dilu- 
tion of 0.001, which corresponds to a SiC column density 
cr = 0.00015 g cm" 2 . From the measured transmission T 
we derived the mass absorption coefficient k = ln(l/T)/(T. 
Though this paper refers to infrared properties, it was nec- 
essary to follow some of the spectra into the visible range. 
Fortunately, KBr is transparent down to a wavelength of 
about 0.2 /jm. Thus, UV/VIS spectra could be measured 
with the same pellets using an Perkin Elmer Lambda 19 
spectrometer. 

3.2. Spectra of high-temperature SiC 

Figure || shows the absorption spectra of four high- 
temperature products in the phonon band range. The 
three fine-grained powders ESK, Lonza UF-15 and AJ-a 
have very similar bands between about 10 and 13.5 /jm 
with the peak at 11.9-12.0 /zm and two shoulders at 
10.7 /j,m and 12.8 /im. The maximum value of the mass ab- 
sorption coefficient is between 15,000 and 20,000 cm 2 g _1 . 
The profiles and absolute values are also very similar to 
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from the visible range to the phonon band region. The 
feature partly overlaps with the phonon band and appears 
in first approximation as a kind of continuous background 
to the latter. 

At a closer view to the literature, one recognizes that 
such a broad feature or background is also present in the 



SiC s pectra measured by other aut hors (e.g. Borghesi et al 
1985| - ff-SiC, [Koike & Shibai 199g| - sample Ibiden). There- 



10 11 12 13 14 15 16 
Wavelength (|im) 



fore, we have to state that this kind of absorption obvi- 
ously is a common spectral characteristics of many SiC 
pyrolysis products. In contrast to this, we did not observe 
the broad absorption in high-temperature SiC samples 
(see spectrum of AJ-a). It certainly is not restricted to 
the /3-polytype as is demonstrated by the a-SiC samples 



Dresden, Plasmachcm, and Ibiden (Koike & Shibai 1998) 



Fig. 5. Infrared band profiles of four high-temperature 
SiC products, solid line: Lonza UF-15, dashed line: ESK, 
dotted line: AJ alpha, dash-dotted line: Duisburg, en- 
hanced by a factor of 10. 



those reported by Fricdcmann et al. (1981)| and Borghesi 



et al. ( L985) 



The extreme broadness of the feature excludes lattice 
vibrations as the carriers of this absorption and strongly 
favours collective vibrations of free charge carriers (plas- 
mons). In small-particle spectra, free charge carrier ab- 
sorption (similar to phonons) causes strong surface modes 
at frequencies smaller than the plasma frequency (Bohren 
& Huffman 1983). These surface modes are very well 



including the shoulders. These have been at- 



tributed by Friedemann et al. (1981) to signatures of the 
transverse and longitudinal optical (volume) modes. How- 
ever, in the light of Sect. 2. we note that at least for the 
shoulder at 10.7 /im this cannot be true because the posi- 
tion of the LO mode is at shorter wavelength. The shoul- 
ders will be discussed further in the following sections. 
Additionally we note here that there is only very weak 
continuum absorption at least at shorter wavelengths. To 
our experience, the continuum is mostly correlated to the 
degree of free carbon pollution in the sample. 

The spectrum of sample Duisburg is characterized by 
a more or less uniform level of relatively small extinction 
with a sharp minimum at 9.6 /im. This behaviour is de- 
termined by scattering and is typical for large grains with 
diameters of 10 fim and more (Andersen et al. 1998). If, as 
indicated by studies of meteoritic SiC, large SiC grains oc- 
cur in the interstellar medium, their different spectral ap- 
pearance would prevent them from infrared detection. For 
the even coarser powders Piesteritz (green) and (black) 
with grain sizes of more than 100 /mi, the measured spec- 
tra are similar to the one of the Duisburg sample but, as 
expected, give even less. Size effects will be treated further 
in Sect. 3.4. 

3.3. Spectra of pyrolysis products 

Figure ^ shows the absorption spectra of the pyrolysis 
products and, for comparison, of one high-temperature 
SiC. Most of the spectra have been measured down to 
a wavelength of 0.5 /Ltm because a large part of them un- 
expectedly revealed an additional strong and extremely 
broad absorption feature, which extends in some cases 



known for metal particles (Kreibig 1996) where they occur 
at visible wavelengths but also for doped semiconductors 
( [Ricdcr et al. 1972| ; |Yamamoto et al. 1985| ). 

The free charge carriers may have been introduced to 
the SiC by impurities, as has been already briefly dis- 
cussed in Sect. 2.4. One possibility is the incorporation 
of nitrogen which, because of its five valence electrons, 
may produce donor levels at different energy depths in the 
band structure ( |Chen et al. 1997 ). For the samples Riga, 
Plasmachem and Dresden the elemental analysis clearly 
detected a nitrogen content. The corresponding spectra 
are characterized by an especially high plasma frequency. 
Since the plasma frequency of a free-electron gas is pro- 
portional to the square root of the charge carrier density, 
this supports the explanation of the broad absorption fea- 
tures by nitrogen incorporation. As a test of the nitrogen 
hypothesis we prepared a nitrogen-doped SiC sample by 
adding a small percentage of NH3 to the reactants. In- 
deed the doped sample (Jena-N) shows the free charge 
carrier absorption in contrast to the undoped one (Jena). 
These investigations will be the subject of a future paper 
(Clement et al., in prep.). 

The phonon band profiles of the pyrolysis products 
show a large variety of shapes, which are obviously not re- 
lated to the polytype. The samples which do not possess 
the free charge carrier feature but only the typical scatter- 
ing increase towards short wavelengths (Jena and Berlin) , 
have phonon bands relatively similar to those of the high- 
temperature powders (compare AJ-a). A more detailed 
look on the band profiles, however, reveals that the pyrol- 
ysis products show the band maximum at 12.2-12.3 /Ltm, 
which is significantly longward from the AJ-a peak. Fur- 
thermore, the typical shoulders of the high-temperature 
SiC are lacking. The Jena-N sample, although it shows 
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the free charge carrier absorption, has a similar phonon 
band feature. The small band at about 8-9 /in seen in 
this and the Plasmachem spectrum is due to surface oxi- 
dation. The Berlin spectrum, additionally, shows a back- 
ground absorption of about 2000 cm 2 g _1 which is due to 
the free carbon excess (see Table ||). One of the most inter- 
esting phonon features is the double-peaked one of the AJ 
whiskers, which will be discussed in the following section. 



The phonon bands of the Riga and Plasmachem sam- 
ples are very weak. The reason for this is totally un- 
clear because the purity of these samples is better than 
the one of, e.g., sample Dresden. One possible explana- 
tion for the weakness of the phonon band might be the 
interaction with the plasmon. This interaction is rather 
obvious for the samples Japan, AJ-/3, and Starck where 
the free charge carrier feature occurs at relatively large 
wavelengths. These spectra show round-shaped phonon 
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Fig. 7. Infrared band profiles of different size fractions of 
the product AJ-a: <0.5 /mi (solid line), 0.5-1 /im (dotted 
line), 1-2 /im (dashed line) and >2 /im (dash-dotted line) 



bands peaking at very short wavelengths (between 11.1 
and 11.4 /im). The sharp depressions at 10.2 and 12.8 /im 
limiting the spectral ranges of the phonon bands, and 
the absorption maxima at wavelengths larger than 13 /im 
point to the theoretically expected coupling of the LO 



phono n to th e plasmon (Sasaki et al. 1989; DiGregorio & 
Furtak l 1993j see also Sect. 2.4). 



Consequently, free charge carriers in SiC grains do not 
only cause near-infrared absorption, but also determine 
the apparent shape of the phonon band if the plasmon res- 
onance is relatively close to the phonon frequency. Both 
facts need to be considered for cosmic SiC grains, espe- 
cially because nearly 1 at.% of nitrogen is contained in 
the meteoritic presolar SiC (Hoppe et al. 1994). 



3.4- Size, shape, and matrix effects 

Size effects on the appearance of the SiC infrared ex- 



tinction have been recently studied by Andersen et al 



(1998) focussing on relatively large grains where scatter- 
ing is dominating over absorption. We focus here on much 
smaller size fractions which we separated by sedimentation 
in acetone from the sample AJ-a with an original maxi- 
mum grain size of 2 /im according to the producer. The 
band profiles of four size fractions are shown in Fig.^, 
they may be compared with the spectrum of the origi- 
nal sample in Fig.[|. The main change in the spectra of 
the different size fractions concerns the long-wavelength 
shoulder at about 12.8 /im which obviously disappears for 
the smallest size fraction. The same effect wc observed 
for the other high-temperature samples. Consequently, we 



can be sure that also the shoulders found by Friedcmann 
et al. (1981) and later authors around the position of the 
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Fig. 8. Comparison of measured infrared band profiles of 
four samples (AJ-a, Berlin, A J whiskers, A 3-/3) with cal- 
culated CDE profiles for /3-SiC taken from Fig. || 



absorption (Bohrcn & Huffman 1983) and will vanish for 



very small particles. 

After correcting for the size effect, the phonon bands 
of our high-temperature samples nicely fit into the LO- 
TO interval 10.2-12.7 /im, but still have the shoulder at 
about 10.7 /xm. Our calculations in Sect. 2.5 do not predict 
such a shoulder at the short-wavelength side of the band. 
In other words, it is not probable that this shoulder is 
caused by the weak anisotropy of the LO-mode frequency 
limiting the absorption band at the short-wavelength side. 
The anisotropy is much stronger for the TO frequencies 
where no shoulder is observed. 

If the optical constants of a material do not explain a 
feature of the band profile, according to Sect. 2 it must 
be caused by the grain morphology. Unfortunately, the- 
ory so far has no exact tool for the calculation of phonon 
band profiles for arbitrary particle shapes. Moreover, from 
the experimental side the situation is complicated by the 
fact that the SiC particles concentrate at the inner grain 
boundaries of the KBr pellet and form agglomerates the 
morphology of which is difficult to study. A simple ap- 



TO frequency indeed are due to contributions from bulk 



proach to shape effects used e.g. by Papoular et al. (1998) 
is the approximation of real grain shapes by ellipsoids. 
In this framework any frequency in a band profile can be 
assigned to a geometrical resonance due to a special ellip- 
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soidal shape characterized by the shape parameter L (see 
Sect. 2.5). The positions of the peak and the shoulder (835 
and 935 cm -1 ) of the smallest size- fraction of Fig. would 
correspond to shape parameters of 0.08 and 0.55, respec- 
tively, which points to more prolate structures. This result 
could mean, that elongated agglomerates dominate in the 
KBr pellet. Considering the limitations of the ellipsoid ap- 
proach, however, such an interpretation remains doubtful. 
We stress here, that for future studies of shape effects it 
is important to separate the influence of agglomeration. 

Figure || compares the already discussed profile of a 
high-temperature product (AJ-a) with some typical band 
profiles of pyrolysis samples as well as with theoretical 
spectra from Sect. 2.5. All the samples belong to the f3- 
SiC polytype. The Berlin spectrum, which is also repre- 
sentative of the Jena samples, still resembles the high- 
temperature spectrum but peaks at a larger wavelength 
(12.2 fim) and has no distinct shoulder. It is the only spec- 
trum which may be related in some way to the theoretical 
CDE spectrum, pointing to the presence of a wide dis- 
tribution of shapes. However, it is obvious that the CDE 
model is insufficient to describe real band profiles. There- 
fore, in the case of SiC it cannot be used e.g. for the deriva- 
tion of optical constants from IR absorption spectra of 
SiC particles. The m-CDE which initially was considered 
to be a more realistic distribution is even worse in repro- 
ducing the measurements, except for the AJ-/3 spectrum. 
However, the band profiles of this and also the Japan 
and Starck samples (comp. Fig.^j), which all peak at very 
short wavelengths (11.1-11.4 /an), obviously are strongly 
influenced by the phonon-plasmon coupling so that a dis- 
cussion in terms of shape effects alone is meaningless. It 
should be noted that the /3-SiC spectrum published by 
Borghesi et al. (1985) , which has been used for fits of as- 
tronomical spectra, seems to belong to the same class. 

The most striking example of shape effects is certainly 
the AJ whisker spectrum. Similar spectra of SiC whiskers 



have a lready been published at different p laces (e.g. Pultz 



& Her f 1966t |DiGregorio fc Furtak 1993|) . All these spec- 
tra are characterized by two peaks, predicted by the theory 
to occur at the TO frequency where the shape parameter 
is zero (796 cm -1 or 12.6 /im for /3-SiC) and at the fre- 
quency where Re(e) =-e m (shape parameter 5, 929 cm -1 
or 10.8 //m in KBr) ( Bohrcn fc Huffman 1983| ). These po- 
sitions are reproduced in our spectrum within 0.2 fim. The 
deviations probably are due to the relatively large size of 
the whiskers. 

All band profiles discussed so far have been measured 
with particles embedded in a KBr matrix which strongly 
influences the band shape. The matrix effects indeed are a 
big problem for comparisons of laboratory to astronomical 
spectra, since at the moment we are not aware of a good 
method to store any given powder sample in vacuum or in 
a gas for laboratory spectroscopy. Since theory is not in 
the position to calculate spectra for arbitrary grain shapes, 
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Fig. 9. Comparison of measurements on particles embed- 
ded in a KBr pellet (dotted lines) and infrared microscope 
measurements on SiC material supported by a NaCl win- 
dow (dashed lines). The solid curves represent the com- 
putational correction of the matrix effect according to the 
method by Papoular et al. (1998) 



it does not provide an exact method for correcting matrix 
effects either. 

Papoular et al. (1998) recently published the already 
mentioned approach, which is valid under the assumptions 
that the particles have a distribution of ellipsoidal shapes, 
and that the band can be described by a Lorentzian os- 
cillator. This method works perfectly in test calculations 
with simulated spectra (if one corrects the spectral inten- 
sities at the calculated new resonance frequencies for the 
changed density of the latter; this step is not mentioned in 
the paper). It was also demonstrated by these authors that 
the method allows estimates of the measurement results 
in different matrices. It is clear that the predictions of this 
theory are not exact because the theory of surface phonons 
of real grain shapes (if a practicable one exists) is not the 
one valid for ellipsoids. Moreover, anisotropy still has not 
been included in the ellipsoid approach. Nevertheless, the 
method is useful for the interpretation of SiC spectra, and 
should definitely be used for approximate matrix correc- 
tion instead of the formerly applied shifting of the whole 
feature. 
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We have compared the predictions of Papoular's 
method with measurements which nearly avoid the ma- 
trix effect. These measurements have been performed by 
means of an infrared microscope with samples dispersed 
on a NaCl substrate (Fig.^). Thus, these grains are not 
fully but mainly surrounded by air. An important disad- 
vantage of this method, which at the moment excludes it 
from being the solution to the matrix problem, is that the 
amount of material in the microscope aperture remains 
unknown. Therefore, the measurement is not quantitative 
but it reveals the shape of the band nearly wit hout a ma- 
trix effect (compare also Andersen et al. 1998 ). 

For three representative samples, Fig. [j] displays both 
the measurements in KBr and with the microscope and, 
additionally, the matrix-corrected spectra in the interval 
lulo^to- For the calculations, the oscillator data of the 
3C polytype and of the 6H polytype perpendicular to the 
c-axis (sample ESK) have been used since the latter dom- 
inate in a hexagonal crystal by 2:1 over those parallel to 
the c-axis. The dependence of the calculated band profile 
on these parameters (within reasonable limits) is not very 
strong. The figure shows that the microscope measure- 
ment confirms the blue-shift of the band center predicted 
by the calculation. This trend is expected because in a 
matrix of smaller refractive index the surface resonances 
have to shift to wavelengths where Re(e) of the particles 
is also smaller. The peak position and the shape of the 
profiles, however, are in only one case reproduced with 
satisfactory accuracy. 

Besides the already discussed inadequacy of the theory, 
certainly experimental reasons contribute to the differ- 
ences. One of them should be the residual influence of the 
substrate, which is known from measurements on metal 



particles (Kreibig 1996). However, using a silicon substrate 
that has a much higher refractive index (Si: n=3.42, NaCl: 
n=1.5), we found no systematic change of the band pro- 
files compared to the results on NaCl substrate. Another 
reason probably is the grain clustering that occurs both 
in the KBr pellets and on the NaCl substrate, but may 
result in different morphologies. 

So at the moment, unfortunately, there is no perfect 
way to solve the matrix problem. Laboratory spectra of 
SiC and other particles with very strong vibrational or 
electronic bands should be compared to astronomical spec- 
tra only very carefully. For SiC we recommend the method 
according to Papoular et al. (1998) at least for a check of 
the possible matrix effect. Hopefully, also for this prob- 
lem future measurements of non-agglomerated grains will 
bring more clarity. 



4. Amorphous SiC 

4--1- Preparation 

Amor phous SiC was obtained by the polymer pyrolysis 
route ( Clement 1996 ). We started from the organometallic 
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Fig. 10. Infrared band profiles of polycarbosilane (PCS) 
heated to different temperatures. From the top: l:room 
temperature, 2:873K, 3:973K, 4: 1123 K, 5: 1273 K. 



precursor dimethylpolysilane (-Si(CH3)2) n (Strem Chem- 
icals) which first was converted to polycarbosilane (PCS ) 
(-CH 2 -SiHCH 3 -CH 2 -Si(CH 3 ) 2 )„ flYajima et al. 1976| ). 
PCS reacts with oxygen and air moisture very quickly to 
silicon dioxide especially at higher temperatures. To pre- 
vent this reaction all preparational steps were performed 
in a glove box and the heating process under dried argon 
atmosphere. 

For the different analyses we have produced the amor- 
phous SiC both in film and powder forms. To obtain the 
films, a thick sheet of PCS was driven on Si substrates. 
For the powders the PCS was kept in an alumina boat. 
The samples were heated slowly in a quartz tube to tem- 
peratures between 873 K and 1323 K which were hold for 
one hour. 

Figure |l0| shows the development of the IR transmis- 
sion spectra of the film samples with temperature. At tem- 
peratures lower than 1023 K the spectra still show fea- 
tures of the organic groups of the precursor PCS (see es- 
pecially the feature at 4.8 ^m). Between 1073 and 1173 K, 
these features have vanished and one broad single vibra- 
tional band from 9.3 to 17.6 fim peaking at 12.65 /xm 
with shoulders at about 12 and 14 /im is observed. X- 
ray diffractometry of the corresponding powder-samples 
revealed only one broad Bragg reflex in the range from 
30 to 40 ° 20 (CuKa radiation). Therefore, these samples 
can be considered to be in the amorphous state. Examina- 
tion by energy-dispersive X-ray analysis gives about the 
stoichiometric ratio of silicon to carbon atoms. 

At temperatures higher than 1223 K the infrared band 
becomes sharper, indicating the beginning crystallization 
of the sample (Fig.|l^). At 1273 K the band ranges from 
10.2 to 17.6 /im with the peak position at the same wave- 
length as the amorphous samples . The X-ray diffrac- 
tograms of these samples show narrower Bragg reflexes 
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and a spike at the position of the (111) reflex of cubic 
SiC. 



4-2. Optical data 

The sample annealed at 1123 K is considered to be the 
best amorphous material and is chosen for the derivation 
of optical constants. The film geometry is ideal for this 
purpose because it can be modelled more exactly than a 
realistic particle shape, for instance. Additional measure- 
ments of the produced powders could not be performed 
because the production of a sufficient volume of submi- 
cron particles failed. 

The layer thickness of the chosen sample, which is 
needed for the derivation of the optical data from the 
IR transmission spectrum, has been determined from 
reflectance spectroscopy in the UV and visible spectral 
ranges. The UV/VIS spectrum shows an interference 
pattern due to multiple reflection in the SiC layer which 
allowes the simultaneous fitting of the dielectric function e 
of the layer material as well as of the layer thickness d. The 
dielectric function was modelled according to Eq. (|l|) by a 
constant £oo = 3.81 and by two Lorentzian oscillators with 
uj to1 =48250 cm- 1 , u) PA = 57819 cm" 1 , 71 = 12935 cm' 1 
and wto,2 = 36185 cm -1 , w^ 2 = 28492 cm -1 , 
72 = 11732 cm -1 , respectively. The film thickness 
and the real DF in the near infrared which is used later 
as the £oo value for the evaluation of the mid-infrared 
transmittance spectrum were determined to be d = 65 
nm and £oo,ir = 5.84, respectively. The latter value is by 
about 10 % smaller than the values for the crystalline 
modifications, which is explained by the usually smaller 
density of amorphous materials. 

For the calculation of the complex refractive index 
from the IR transmission spectrum of the thin film on sil- 
icon substrate (1123K), a Lorentz-oscillator fit has been 
used, too. For an exact fit of the band range we applied 
a large number (37) of oscillators at different wavelength 
positions. The need for so many oscillators is a strong con- 
trast to the properties of " crystalline spectra" and reflects 
the presence of a wide distribution of vibrational modes 
due to the different bonding lengths in the amorphous 
material (comp. Sect. 2.4). The resulting complex refrac- 
tive index is shown in Fig.jll]. Compared to the bands of 
crystalline materials it is, as expected, much broader. The 
same holds for the small-particle absorption spectra calcu- 



lated f fom these optical constants (Fig." 12, comp. Figs. 5, 
||, 0). The accompaniing "features" at 7 and 22 fim are 
artefacts caused by the limited accuracy of the substrate 



compe: isation. 



5. Comparison with observations 

The goal of this paper is to give a detailed laboratory basis 
for an interpretation of astronomical SiC profiles. We do 
not aim at analyzing these profiles here, this should be the 
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Fig. 11. The complex refractive index (real part — dotted 
line, imaginary part = solid line) calculated from the amor- 
phous SiC layer hold at 1123 K. 
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Fig. 12. The small-particle absorption spectra calculated 
from the optical constants (Fig.|ll|): for spherical grain 
shape (solid line) , for a continous distribution of ellipsoidal 
grain shapes (dotted). 

task after more high-resolution ISO spectra have become 
available. However, we will discuss the relevance of some 
of our results for astronomical infrared spectroscopy. 



If the interpretation of IRAS spectra by Goebel et al 



(1995) and Baron ct al. (1987) in terms of an a:C-H con 



tribution is correct, then only the stars with relatively 



high continuum temperatures (class 4 in Goebel et al 



1995) have features of purely SiC origin. These features 



peak always around 11.3 /im. This is also supported by 



the UKIRT observations by 3peck et al. (1997), which are 



probably the pre-ISO carbon star spectra of the best qual- 
ity so far. In these data, for continuum temperatures above 
1000 K (e.g. of R For, IRC+50096, UU Aur, V Cyg) fea- 
tures dominate, which peak around 11.3 /xm with a nearly 
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symmetric profile and a width of around 1.75 /im. The 
profiles are not entirely symmetric, because there is an 
indication of a shoulder at about 11.8 fim. 

In contrast to this, most of our laboratory spectra 
show band maxima at relatively long wavelengths, e.g. the 
high-temperature products at 11.9 /im and most pyroly- 
sis products around 12.2 /im. Whether the matrix effect 
could have transformed band profiles peaking originally 
at 11.3 [im into the measured ones is still open, according 
to our infrared microscope measurements in Sect. 3.4 we 
rather doubt it. It is much more probable that morpholog- 
ical effects of the grains in the KBr pellets determine the 
measured laboratory profiles. However, up to now this has 
only been shown for a few samples with extreme geome- 
tries (whiskers), since at the moment the agglomeration 
state of the laboratory samples cannot be taken into ac- 
count. 

The only laboratory samples which show the phonon 
features at about the correct position (11.1-11.4 /im) and 
with the correct width are those for which the bands are 
influenced by the phonon-plasmon coupling. To decide if 
the circumstellar SiC is of this type, there is in principle 
the possibility to search for an indication of the very broad 
absorption at adjacent wavelengths or for sharp minima 
limiting the phonon band region. 

Whether such an identification will be possible in the 
future, depends also on the knowledge about other fea- 
tures of the astronomical spectra, especially those caused 
by a-C:H and molecular species. Their contributions have 
to be known for an exact determination of the base line 
of the SiC band. Another problem of the identification 
may be self-absorption, which has already been addressed 
by |5peck et al. (199T) who used it for fitting their band 
profiles. Self-absorption may change the band profile dras- 
tically and produce relatively broad profiles with flat tops 
as well as maxima at longer wavelengths. This might 
also be another possible explanation for the peak shift to 
11.9 /im observed for the lower continuum temperatures, 
since cooler evelopes should have larger optical depths. 
Indeed, these observed features are relatively broad which 
may be difficult to explain with our data. However, they 
still are considerably narrower than the amorphous SiC 
feature we have measured (Sect. 4). Therefore, amorphous 
SiC can be ruled out as an analog of circumstellar SiC. 

Another observational fact of importance might be the 
mentioned shoulder at about 11.8 /urn. If this is not due 
to a:C-H absorption, it could be indicative of larger grain 
sizes of a few micron (comp. Sect. 3.4). In this light, the 
whole observed feature shift might also reflect a growing 
grain size which would be understandable since lower con- 
tinuum temperatures are correlated with higher mass-loss 
rates ( Skinner fc Whitmore 1988 ). The observed weaken- 
ing of the feature is also expected from a growth in grain 
size. 



6. Conclusions 

In our discussion of the optical properties of various SiC 
crystal types in Sect. 2 we have shown that the optical 
constants characterizing the phonon excitation by infrared 
radiation are too similar to cause significant spectral dif- 
ferences in small particle spectra. Nevertheless, in Sect. 3 
we have presented spectra of considerable variety which 
is definitely not correlated to the polytype. Consequently, 
morphology and impurities of the material are the main 
factors determining the band profile, although there may 
still be indirect effects of the polytype, e.g. via the grain 
shape. This finding opens a wide field for future studies 
which are necessary for the understanding of the cosmic 
silicon carbide dust. 

The points which have already been clarified or ad- 
dressed in this paper are: 

1. The long- wavelength shoulder of the SiC bands mea- 
sured by Friedemann et al. (1981) is a size effect. 

2. The short- wavelength shoulder seen in these spectra is 
not positioned at ujlo- It may be caused by a shape 
effect but this is unclear so far. 

3. Free charge carriers strongly influence the profile of the 
phonon band via coupling to the LO phonon. 

4. The free charge carriers themselves are a source of 
strong absorption in the mid-IR. They are very com- 
mon among SiC powders. Probably, most of the free 
charge carriers are produced by impurities. At least 
one of the possible dopants is nitrogen which is also 
found in the presolar meteoritic SiC grains. 

5. The crystal structure (polytype) is not a strong factor 
for the band profile, at least not in the direct way via 
the dielectric function. 

6. The infrared band profile of amorphous SiC is too wide 
to fit the 11+ feature in circumstellar envelopes of car- 
bon stars correctly. 
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